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We discuss cosmological perturbation theory at third order, deriving the gauge transformation 
rules for metric and matter perturbations, and constructing third order gauge invariant quantities. 
We present the Einstein tensor components, the evolution equations for a perfect fluid, and the 
Klein-Gordon equation at third order, including scalar, vector and tensor perturbations. In doing 
' so, we also give all second order tensor components and evolution equations in full generality. 
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I. INTRODUCTION 



Anisotropies in the cosmic microwave background and large scale structure provide compelling evidence that the 
universe is not truly homogeneous and isotropic and therefore cannot be described completely by the Friedmann- 
Robertson- Walker (FRW) spacctime. Since general relativity is highly nonlinear, it is extremely difficult to find 
exact solutions which would allow for all the inhomogeneities and anisotropies present in the real universe. Thus one 
v««/ ■ usually resorts to perturbative techniques, considering a homogeneous FRW background and adding inhomogeneous 
' perturbations. This is the basis of cosmological perturbation theory, which has become a cornerstone of modern 
P ,! cosmology in the last half century (sec for example Refs. [J, [3, i, i, [I B, 0, B B IHI, [13, El ) ■ 
' ' Linear perturbation theory is now an essential tool in the cosmologist's toolbox, and extending the theory beyond 
^ , linear order is at present a rapidly expanding area of research: recent years have seen second order theory gaining the 

■ attention it duly deserves and becoming more mature. However, one does not have to stop there: it is worthwhile and 
Ci , feasible to consider perturbation theory even beyond second order, and in this paper we study aspects of cosmological 

perturbation theory at third order. 

Each order in perturbation theory reveals different though complementary aspects of the underlying fully non-linear 
^ . theory, and also gives us access to information of different data sets or to extract more information from the same data 

■ set. First order theory allows us to model the large scale structure of the universe and to calculate the distribution of 
Tij" , anisotropies in the CMB. Second order theory then enables us to access even more information from these data sets, 
0^ ' by calculating higher order observables such as the bispectrum. Third order reveals yet even more information, by 

i calculating, for example, the trispectrum. 
^\ • Whereas at linear order the different types of perturbations, denoted by their transformation behaviour on 

■ three-hypersurfaces as scalar, vector and tensor, decouple, this is no longer the case beyond first order. At second 
0^ . order we find, for example in the energy conservation equation, couplings between first order scalar and vector 

perturbations through gradients, and also the coupling of first order tensor perturbations to each other. At third 
order a new coupling occurs in the energy conservation equation, namely the coupling of scalar perturbations to 
tensor perturbations. This will allow for the calculation of yet another different observational signature, highlighting 
another aspect of the underlying full theory. 



There has already been some work on third order theory. For example, Ref. |14l . Il5| considered third order 
perturbations of pressureless irrotational fluids as "pure" general relativistic correction terms to second order 
quantities. The calculations focused on the temporal comoving gauge, allowing the authors to consider only second 
order geometric and energy-momentum components, and neglected vector perturbations. Ref. [iGj includes a study 
of third order perturbations with application to the trispectrum in the two-field ckpyrotic scenario in the large scale 
limit. There has also been reference in the literature of the need to extend perturbation theory beyond second order. 
For example, in Ref. [l3| UV divergences in the Raychaudhuri equation are found when considering backreaction from 
averaging perturbations to second order. The authors state that these divergences may be removed by extending 
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perturbation theory to third, or higher, orders. 

In this paper, we develop the essential tools for third order perturbation theory, such as the gauge transformation 
rules for different types of perturbation, and construct gauge invariant quantities at third order. We consider perfect 
fluids with non-zero pressure, including all types of perturbation, namely, scalar, vector and tensor perturbations. In 
particular allowing for vector perturbations is crucial for realistic higher order studies, as we have recently shown that 
vorticity is generated at second order in the perturbations in all models employing non-barotropic fluids p^ . Hence 
studying irrotational fluids at higher order will only give partial insight into the underlying physics. We present 
the energy and momentum conservation equations for such a fluid, and also give the components of the perturbed 
Einstein tensor, up to third order. All equations are given without fixing a gauge. We also give the Klein Gordon 
equation for the scalar field minimally coupled to gravity at third order in cosmological perturbation theory. 

The paper is organised as follows: in the next section we give our definitions, followed by the gauge transformation 
rules in section IIIII We construct gauge invariant variables at third order in section IIV| present the governing 
evolution and constraint equations in section |V] and conclude with a discussion in section IVII 

In this paper, we use conformal time, 77, throughout, denoting derivatives with respect to conformal time with a 
prime. The scale factor is a, the comoving Hubble parameter is 7i = a'/a and the physical Hubble rate is H = a' /a^. 
Greek indices, fJ-,v,X run from 0, ...,3, while lower case Latin indices, i,j,k, take the value 1,2, or 3. Covariant 
derivatives are denoted by a semi-colon, partial derivatives by a comma. The order of the perturbations is usually 
denoted with a subscript immediately after a perturbed quantity, when it is split up order by order. We assume a 
spatially flat FRW background throughout. 



II. DEFINITIONS 

The components of the covariant metric tensor are 

500 = -a^(l + 2</'), goi^a^B,, = a2(% 2Cy) , (2.1) 

where we take the three dimensional background space to be flat with metric Sij, and include scalar, vector and 
tensor perturbations. Here, (f> is the lapse function and Bi is the vector perturbation which can be split into the 
gradient of a scalar, B^i, and a divergence-free vector. Si} Cij is the spatial part of the metric perturbation, which 
can be decomposed as Cij = —ipSij + E^ij + F(^i^j) -\- \hij. Note that in the uniform curvature gauge, this reduces to 
Cij = \hij, where hij is the tensor, or gravitational wave, perturbation. 

The contravariant metric components are obtained by imposing the constraint g^i^g'^^ = 6^^, to the appropriate 
order. To third order this gives 

.g"° = -4 (1 - 20 + - 803 - BkB'' + ^BkB'' + 2B'B^Ci,) , (2.2) 



g°' = ^ (S* - 2(f)B' - 2SfeC'=' + 4(f>^B' + ABkC'^'cj) + AC'^^Cj'Bk - B'^BkB') , (2.3) 
gij ^ J_ tgU _ + AC'^Ck^ - B'B^ + 2(j)B'B^ - ^C'^C'^^Cki + 2B'C''^Bk + 2BkB'C'^) . (2.4) 
The perturbations can be split order by order as, e.g.. 



(2.5) 



but we do not explicitly perform such a split here in order to keep expressions compact. For example, expanding the 
— component fully order by order gives 

= (1 - 201 - 02 - ^03 + 0? + 80102 - 80? - BifcSi^ - B2fc-Bi^ + 40iBi^Sifc + 2Bi'B^^Cu^ , (2.6) 



^ Note that our sign convention differs to that in Ref . |15|| , by a minus sign in front of Bi . 
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(2.8) 



which when compared to Eq. (|2.2p illustrates why we do not split perturbations where possible in this paper. 

The energy momentum tensor for a perfect fluid is 

T^, = (p + P)u^u, + PS^', , (2.7) 
where P is the pressure, p is the energy density and u'^ is the fluid four velocity, defined as 

The fluid velocity is subject to the constraint 

u''Uf, = -l, (2.9) 
and, to third order in the perturbations, has components 

u' = -v\ (2.10) 
a 

u" = i (^1 - + - ^cf,^ + ^Vkv'' + WfeS^- + Ck,v''iP - 2<j>v''Bk - , (2.11) 

u^^a U, + B,- 4>B, + 2C,kV^ + ^B,^^ + ]^B,v''vk + B.v^bA , (2.12) 



Wo = -a 1 + - -0-^ + -0^ + -w'^Wfc + (tyvkv" + (7^^!;%^ . (2.13) 



III. GAUGE TRANSFORMATIONS 

We use here the active approach to deal with gauge transformations, where the exponential map is the starting 
point (we refer the reader to the recent reviews on perturbation theory [13, |Tl| for detailed discussions of gauge issues 
and related matters in this section ). This allows us to immediately write down how a tensor T transforms once the 
generator of the gauge transformation, has been specified. The exponential map is then 

T = e^"«T, (3.1) 

where £^ denotes the Lie derivative with respect to The vector field generating the transformation, , is up to 
third order 

^ + + \^^e, + 0[e^) . (3.2) 

The exponential map can be readily expanded and, taking into account Eq. (|3.2[) . gives 

exp(i:e) = 1 + e£^, + ]^e^£l + ]^e'£^, + ^e^£^, + £l + -/£^,£i, + ^e^% + ■ ■ ■ (3-3) 

where we have kept terms up to 0(e^). Splitting the tensor T up to third order, as given in Eq. (|2.5|) . and collecting 
terms of like order in e we find that tensorial quantities transform at zcroth, first, second, and third order, respectively, 
as [il[ll 

To = To, (3.4) 
e5Ti = e^Ti +ei;5,To, (3.5) 
£2^2 - (^T2 + £^,To + 4^ To + 2%,5Ti) , (3.6) 

e'<5T3 = e3|5T3+ + + + ^%%^To + 3(4^+%)5Ti+3%(5T2|. (3.7) 

Applying the map (j3.ip to the coordinate functions a;^, we get the following relationship between the coordinates 
of a point p and a point q 

x^'{q) = e^^^\^ x^'{p), (3.8) 
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where we have used the fact that when acting on scalars £j = ^^{d/dx^) and the partial derivatives are evaluated at 
p. The left-hand-side and the right-hand-side of Eq. p.Sp are evaluated at different points. Equation (|3.8p can then 
be expanded up to third order as 

x^{q) = x^(p) + eer(p) + ie2(CbKi'^b)+e2^(p)) 

+ (Ca (P)^'(P) + e,,Mi2ip)) , (3.9) 
relating the coordinates of the points p and q. 

A. Four-scalar s 

We now turn to four-scalars, using here the energy density p as an example, and study its behaviour under a gauge 
transformation up to third order. It can be readily expanded as 

p = po + ^Pi + ^5p2 + ^<5p3 , (3.10) 
and we shall now detail the transformation order by order. 

1. First order 

The Lie derivative of a scalar is simply given by 

£^p = ep,x. (3.11) 

Before we can study the transformation behaviour of the perturbations at first order, we split the generating vector 
into a scalar temporal part ai and a spatial scalar and vector part, /?i and , according to 

= (ai,/3i;+7iO , (3.12) 

where the vector part is divergence-free {dk^^ = 0). 

Under a first-order transformation a four scalar, here the energy density, p, then transforms, from Eqs. (j3.5[) and 
(laHl) as, 

5pi^5pi+ p'^ai. (3.13) 

The first-order density perturbation can be made gauge-invariant by prescribing the first order temporal gauge or 
time slicing, ai. 

2. Second order 

At second order, as at first order, we split the generating vector ^2 i^^to ^ scalar temporal and scalar and vector 
spatial part, as 

e2={o^2.p2:+l2) ^ (3.14) 

where the vector part is divergence-free {d^'-^^ = 0). We then find from Eqs. (|3.6p and (|3.1ip that a four scalar 
transforms as 

5p2 = 5p2 + PQa2 + ai (/Oo'ai + Po"i' + 2(^/5i') + (2(5pi + Poai) (/Ji/" -I- 7i''') . (3.15) 

We now need to specify the time slicing at first and second order, and also the spatial gauge or threading at first 
order, in order to render the second order density perturbation gauge- invariant. 
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3. Third order 



At third order we also split the generating vector ^3 into a scalar temporal and scalar and vector spatial part, as 

= («2,/32,' + 730 , (3-16) 

where the vector part is again divergence-free {dkj^^ = 0). Wc then find from Eqs. (|3.7p and (|3.1ip that the energy 
density transforms as 

Sp3 = + (^£^^ + £^^^ + ^£^^£^2 + PQ + 3 + £^.,) Spi + 3£^^Sp2 , (3.17) 

which gives 

Sp3 = 6p3 + p'^as + Po'al + Sp'^aiai^x^i + Po (ai,A/3Ci + "i.A^i, /?) 

+3p^'aia2 + Pol {a2Mi + aLA^a) + 3 {Spi^xp^^ + Sp^^^^l ^) + SSp^x^^ + 3Sp2.x^i • (3.18) 

Similar to the second order case, we need to specify the time slicings (at all orders), and also the spatial gauge or 
threading at first and second order, in order to render the third order density perturbation gauge-invariant. 



B. The metric tensor 



We now give the transformation behaviour of the metric tensor. As above, we refer to Refs. [T3| and [Tl| for details 
on the first and second order calculations, and focus the detailed calculation on the third order derivation. The 
starting point is again the Lie derivative, which for a covariant tensor is given by 

£^9^.. = Sm-.a^^ + 5MAe.\ + .9A.e.\ • (3.19) 



1. First and second order 



The metric tensor transforms at first order, from Eqs. (|3.6p and (|3.19p . as 



Off^i/ +5^i.,A?l +5pA?l ,u+9xu^l , 



(3.20) 



and at second order as 



AO) 



JO) 



.(0) 



„(0)^A pa 



+f ^lA (^1 ,iya^l + ^1 ,Q^" 1^) + 9\J (^1 ,fj.a^l + ,a^l, 



(0) 



(3.21) 



Using the method detailed in Refs. [T3] and [Tl| we can extract from the general expressions Eqs. (|3.20p and (|3.2ip 
the transformations behaviour of a particular metric function. We focus here on the curvature perturbation ip, and 
find that it transforms at first order as [1, 0] 



4'i — Hai 



(3.22) 



At second order we get after some calculation 



V'2 = i'2-na2 - -x-^k 



(3.24) 



^ The general expression given above in Eq. I)3.21|l gives the transformation for C2ij , namely 

2C2,j = 2C2rj + 2na2hj + + + ■^ij , (3.23) 
and hence intermediate steps are necessary to extract the transformation of a particular component (see Refs. lid and llll for details). 
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where X2ij contains terms quadratic in the first order perturbations. For ease of presentation, on taking only scalar 
perturbations into account and working on large scales, where we can neglect gradients, X2ij takes the simple form 



2ij 



and Eq. (f?^ reduces to 



V'2 = "02 - "^012 



a 
a 



al - Haia[ + 2ai (V-i - Hipi) 



(3.25) 



(3.26) 



2. Third order 

As above in the case of the transformation behaviour of a four-scalar at third order, the change under a gauge 
transformation of a two-tensor can be found applying the same methods as at second order. We therefore find that 
the metric tensor transforms at third order, from Eqs. (|3.7p and (|3.19p . as 



3 ,/i 



+3 



(2) fA 



.(2)fA 



+ \ 9fj.\,i3^2,v + 9x^,(3^2,11 + 5pA ^2,/3i/ + .9ai/ ?2,/3a' J ^ 1 + \^.9/3i/,A?2 + 5ai/ ?2,/3 J 
_L tfA , „(0)^A (JO) >A , (0) .A , (0).A , (0).A \ 

+ \^5^i/,ASl + 5ai/ S1,/3 J ?2,zy + l^ff/3iy,ASl + ffAi. ?l,/3 J ^2,a» 

, J„(0) COf'S^A , o„(0) tXfrafrP , (0) ,A ^a;j/3 , ^A >a 

+ S 9fj.^,af3X^l ^1 ^1 + •^9^iv,Xp^l 5l ^l,a + S^iy.A'Sl^Q/J^ 1 ?! + 5^iy,ASl,/35 

+'^3m/3,A (Cl Sf^aC^.iy + C?Sl Sf,Qiy) + 3g^a,ASf Sl,/?^?,!^ + ^S^J./JA^f Si S" i/ 

tafPtX \ fA I fOif^' fA I oe/^fA /-Q I f^'fA fa 

SI ^1 ^l,al3v + S1,qS1,/3?1,i/ + ?1 ^l,a^l,l3y + ^Sl Sl,a/3?l,i/ + ?1 ^I.qSI,/?!^ 

Si Sf.aS",! + SrSi Sf,Q^) + ^g^Q.ASf Si'.^Si',^ + 3g!,Q,/3ASf SiS"^ 



fOflSfX I <:P <rA <^q i C^C^ cA i Ot^fA \ f>^ fa 

SI SI Sl.Q/3^ + Sl,aSl.;3Sl,^j + Si Si,qS1,/3^i + ^Sl Sl,a/3Sl,p + Si S1,qS1, 



/3m 



+"ffa/3,ASlSl,,y + JffaA 



Sl,/i (^S?,^Sf,jy + S"iySf,^) + Si (S"/jSl,/3iy + S"i/Sl,/3p) 



(3.27) 



However, in this case it becomes even more obvious than in section [IT] that the expressions at third order are of not 
inconsiderable size. This will also be clear from the Einstein tensor components and the evolution equations given 
below in Section IVl 

Now, following along the same lines as at second order, Eq. p.27p gives the transformation for the spatial part of 
the metric at third order. 



2C3y = 2C3ij + 2T-ia^6ij + 2^3(.jj-) -I- X^ij , 



(3.28) 
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where X^ij contains terms cubic in the first order perturbations. Extracting the curvature perturbation gives 



V>3 = ^3 - -Has - + ^-^X^^^, 



(3.29) 



This expression is general, including scalar, vector, and tensor perturbations and is valid on all scales. However, as 
above we shall detail here only the expression valid for scalar perturbations and large scales, at third order an even 
bigger blessing than at second order, and find that X^ij takes then the simple form 



X- 



3ij = '^a^Sij { — 5 



m- 



2 



H {aia'.2 + a[a2) + 2 h ] aia2 



Hence we finally get for the transformation of 



V'3 



a 



.a 



a 



a 



i'3 + na3+ [— + 3H— I + 3 ( — + ala[ 



H {pL\OL2 + ol-iOfi) + 2 



(3.30) 



(3.31) 



IV. GAUGE-INVARIANT QUANTITIES 



In the previous section we have described how perturbations transform under a gauge shift. We can now use 
these results to construct gauge-invariant quantities, in particular the curvature perturbation on uniform density 
hypersurfaces, C- In this section, as before, we consider only scalar perturbations, and restrict ourselves to the large 
scale limit. 



A. Defining hypersurfaces 



From Eqs. (j3.13|) . (|3.15p . and (j3.18p we find the time slicing defining uniform density hypersurfaces at first, second 
and third order in the large scale limit as 



0L2&P 



Sfh 
P'o ' 

Sp2 , Spi6p[ 
p'o 



I 2 

Po 

1 



'^P'o 



3 {6pi6p'2 + Sp[5p2 



p'o 



4Vr^^ 
Po 



P'o^P'l 



5pi 
p'o 



(4.1) 
(4.2) 

(4.3) 



Similarly, the temporal gauge transformation on uniform curvature hypersurfaces is defined by evaluating Eqs. (j3.22p . 
(|3.26p . and (|3.3ip and gives, at first, second and third order. 



"ifiat — ^ , 



a2flat 



-03 

n 



n H TP ' 
1 



2W 



3-01 ^^2 H 7^ 1- oHipitp2 



H 



H 



ipftp[ fa" 37 
7^4 



(4.4) 
(4.5) 
(4.6) 
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B. Constructing gauge-invariant quantities 

We can now combine the results found so far to get gauge- invariant quantities, and as before choose the density 
perturbation on uniform curvature hypersurfaces and the curvature perturbation on uniform density sHces as examples. 

1. Energy density perturbation on uniform curvature hypersurfaces 

A gauge invariant matter quantity of interest is the perturbation to the energy density on uniform curvature 
hypersurfaces. This is obtained by substituting the temporal gauge transformation components in the uniform density 
gauge (given in section HV Ap into the appropriate transformation equation (that is, Eqs. (j3.13p . (|3.15p or (|3.18p ). 
This gives, at first, second and third order, respectively. 



dpiflat dpi + Po^ , 



(4.7) 



H n 



(4.8) 



Opsflat = 0P3 + Po^ 



"0? r „/ „ / a" \ , / a" 

• 3H) +3p'oH(3— -n 



"' o 'ft 

Po -3po — 
a 



a 



Po 



2p'on-p'o- 



a \ a 



(4.9) 



2. Curvature perturbation on uniform density hypersurfaces 



The curvature perturbation on uniform density hypersurfaces, ^, is defined as 



(4.10) 



This is obtained by substituting the temporal gauge transformation components in the uniform curvature gauge (given 
in section HV A|) into the appropriate transformation equation (that is, Eqs. (|3.22p . ()3.26p or (|3.3ip ). Evaluating this 
on spatially flat hypersurfaces then gives, to first, second and third order, respectively. 



Po 



(4.11) 



-—6p2 
Po 



2H 



SpiSpi' 



P'i 



n 



Po 
P'o 



(4.12) 



Sp3 

p'o 



p'o' 



^^{Sp2'Spi + 6pi'Sp2) 
Po 



Po 



Sp,'Sp^'\2(- + 



—Sp2Spi 
Po 

Spi"" 



Wo 
p'o 



p'i 



-n 



El 
p'o 



m 



Spi'^Spi" 



6^ A '2 a 

-^6pi dpi 



a 



p'^y a 



(4.13) 



There are different definitions of the curvature perturbation present in the literature, depending on different decom- 
positions of the spatial part of the metric tensor. A different definition to the one above, as discussed in e.g. Ref. [lH, 
was used by Maldaccna in Ref. [2l| to calculate the non-gaussianity from single field inflation, and was introduced by 
Salopek and Bond in Ref. [13]. They define the local scale factor a = e", then 



a\i^)^^ = a\ij){l + 2CsB + 2Cs^B + xCsb) ■ 



(4.14) 



^ The sign is taken so as to coincide with f as defined in Ref. [20II . 
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Comparing to the expansion from perturbation theory 



a'(^)(l + 2C), 



(4.15) 



one can obtain the relationship 



C-CSB 



+ CsB + qCsB ■ 



(4.16) 



Splitting this up order by order gives, at second order 



C2SB 



C2-2(Ci)^ 



(4.17) 



and; at third order, 



Note that it is this definition, (|4.18p . of the curvature perturbation 
pre-factors since their perturbative expansion is defined differently'*. 




V. GOVERNING EQUATIONS 



Having constructed gauge invariant quantities up to third order in the previous section, we now turn to the evolution 
and the field equations. 



In this subsection, we give the energy momentum conservation equations for a fluid with non-zero pressure and in 
the presence of scalar, vector and tensor perturbations. The latter generalisation is important since at orders above 
linear order, all types of perturbation are coupled. It is also important to consider rotational fluids, since second 
order vorticity in the early universe is likely non-zero (see, for example, Ref. [3]). This agrees with Ref. p^ . but is 
more general. 

The Bianchi identities imply, through the Einstein field equations which relate the geometry of the space-time to 
its matter content, energy momentum conservation, 



where here denotes the covariant derivative. Substituting the definition for the energy momentum tensor (j2.7p 
into Eq. (jS.ip gives energy conservation (the zero equation) 



A. Fluid conservation equation 



V^T^, = 0, 



(5.1) 



+ (po + Po) [{B' + 2v'){v'^ + B[) + v\^4> - 2C[fi'' + v\C' + 20,,) + mv\B^ + 2vi) 
+ ((5p + bP) [(B' + 2i>*)(w^ + B[) - 2C[f'^ ^ v\C\^, ^- 20,,; + m(v, + B,)) + w\,0 
+ {6p' + 5P'y{B, + Vi) + {5p + 5P),,4>v' + (po + Po)(2C"'^w.Wj - B.v'cf) + v'v,(f>) 




(5.2) 



* It is worth mentioning that fsBi the variable first introduced by Salopek and Bond and then employed for non-gaussianity calculations 
by Maldacena, is extremely Gaussian after slow-roll inflation, as opposed to other f variables which exhibit non-gaussianity, as can be 
seen e.g. from Eq. 1 14. 171 1. We thank the anonymous referee for highlighting this issue. 
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and momentum conservation (the z— component) 

+ {6p + + AU{v, + Bi)) - + P^) \{2B, + Ui)^ - 2C,jV^ 



+ (2CywJ')' - (t)ivl + 2B[ + 20,, + mv,) + (pj, + P(J) + Bjw, + 2B,Bj + 



2Ci 



+ i^v, + AB,)cj)^ 



3„y 



2ni8B, + 3v,) + ^ 



+ (Po + Po) 

+ {v, + 2B,)(v'^B^ - C^'^cf) + {v, + - + i'^' {25,(5; + v'^) 

+ Bj(2B^ + 20,, + v[) + 2U{iH + 2B,){2Bj + v^) + 2Cy (C'^'^ + w^fc - 4H0) 

+ C\,j(B,; + W,) + V,(B; + 0,j) + Vj{B', + 0,,) + v'=(2C,;fej - Cjfe,,) 

+ (Bj,, - B,.j - 2C^)0 + 2akv\j} + ^iv^VJv'y - 2C,jW^'0 + B,(40V - Wj/'c 



= 0. 



(5.3) 



Equation (|5.2p highlights the couphng between tensor and scalar perturbations which occurs only at third order 
(and higher) in perturbation theory. At both linear and second order, no such coupling exists, since the only terms 
coupling the spatial metric perturbation, C'ij, to scalar perturbations contain cither the trace or the divergence of 
Cij and the tensor perturbation, hij is, by definition, transverse and trace-free. However, at third order, terms 
like SpC'^jC^^ occur in the energy conservation equation which, on splitting up order by order and decomposing 
becomes Spih'^^^hi^ . It is clear that this term only shows up at third order and beyond. Thus, as mentioned earlier, 
third order is the lowest order at which all the different types of perturbations couple to one another in the evolution 
equations, which will produce another physical signature of the full theory. 



1. Scalar s 



It will be useful to have energy and momentum conservation equations for only scalar perturbations. These equations 
are obtained by making the appropriate substitutions Cij ~ —Sijip + Ejj, Vi = u,^, and Bi = Bi into the above 
expressions. On doing so, we obtain the energy conservation equation 

Sp' + m{5p + 5P) + (po + ^o)(V2u + V^E - 3^/'') + {5p + 6P){V^v + V^E - Sip') 
+ {dp + 6P),,v; + (po + Pq) [{B' + 2v;){v'.^ + B'J + V^v(t> - 2{i''{i'4^ - V^E) 
~ ipV^E' + E[^jE;^^ + v;(20 + V^E^, - 3^-,,) + 4Hw;(B , + 2w,,) 
+ (dp + SP) [(S; + 2v;){v[^ + B'i) - 2[iP'{i^ ~ V^E) - i^V^E' + E[^^E,''^ 
+ z;;(20,, - 3^,, + V^E + m{v^, + B ,)) + V^w^l + [dp' + 5P')v.'{B^, + w,,) 
+ {5p + SP)^,(l)v^' + [po + Fo)(w>,»'/> - S V - eipv^^vj + E^'^v^,vj) 

- {po + Po) [2E;^v^''E,,jk - 2v>j + B ,5;^. - 2E/E^kj) - 2\/^Eij,kV,'' 

- 2il;{v!'{y'^E^k - 3V'.fc) - 2w^u/' + B'^B,^ + GV-^/'') + ^\7^v{(j)^ - v/vj) 

+ {v; + V^E' + B^^){2B',(j) + 3^/v^,) - t>/'|B;'S,y + u>,y - 20'Bj - S^j'v^j 
+ 0(wj + 2v[^ + SHv^j - 2(t)j ~ + V^E.j) + E[^^v;} 

+ B,^ (B,j- V'' + E[^jv; + B^j(P' + v[j(j))] = , (5.4) 



and the momentum conservation equation 
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{po + Pn){t\^ + S ,)] + {po + Po){<p,^ + m{v,^ + S ,)) + (5F + [{5 p + SP){v,, + B,,)] 
+ [Sp + 5P){^,, + 4H(«,, + B,,)) - (Po + ^o) [(2B . + + "^^v,, - 2B,yZ;/' 
+ (po + Po) [{v,^ + B^,){V\ - SV-' + V'E') - B4<j)' + 8H<j)) + t'/' (v^,, - 8H{^S^j ~ E^,,)) 

- 2(V'w,. + E^.jv,^)' - (j){v[^ + 2B;, + 20,, + AHv^,)] + (pj, + P^) [w/'(B,,t;j + BjU,,; + 2B,,Bj 

+ + 20(^<5y - S,,,)) + {^v,, + 4B,,)0'] + (po + Po) + 25,,) (i;^^/' - (St^' - V^E') 

+ t/jW^E' - E/'^E'^k^ + w,^|2B,,(B^ + v[^) + B,j-(2B;, + 20,, + i/,) + 27^(w,, + 2B,,)(2S,j- + z;,, 

- (30,, - V2£;,,)(B., + w,,) - 2(0>5„ - i;,,j)(V2w + \/^E' - 30;' - 4H0) + z;,,(B;j + 0,,) 
+ vj{B[^ + 0,,) - 20.jw,j + i^,,fejw.'' + V^,jw.j + 2(0'(5,j - B.,j))0 - 20w,,j + 2i;,,fcU,'=j| 



0. 



(5.5) 



Note that, in the above, denotes the Laplacian; i.e. = djd^ . 

Considering the large scale limit, in which spatial gradients vanish, the energy conservation equation becomes 



(5.6) 



Sp' + m{5p + 5P) - 30'(po + Po) - 30'(<5p + 5P) 
- e^b^'ipo + Po) - 600' {Sp + SP) + 1202^' (po + Po) = . 
Splitting up perturbations order by order, this becomes 

Sps' + 3n{Sp3 + SP3) - 30^(po + Po) - 90i((5pi + SPi) - 9i'[i6p2 + SP2) 

- 18(po + Po){MA + 0102) + 7202^; (po + Po) = . (5.7) 
In the uniform curvature gauge, where = 0, this is 

Sp'sn,, + 37^((5p3flat + SPsti^t) = , (5.8) 
and in the uniform density gauge, where Sp — 0, 

mSP^Sp + 3C^(po + Po) + X'^SPiSp + %[SP2Sp - 18(po + Po)(C2Ci + C1C2) - T^ClCiiPo + ^o) = , (5.9) 

with C as defined in section flVB 21 This can be recast in the more familiar form by introducing the (gauge invariant) 
non-adiabatic pressure perturbation. At linear order the pressure perturbation can be expanded as 



dP dP 

SPl = "^^-^l + 'g^^P^ ^ + <?JP\ ■ 



This can be extended to second order [23| and higher by simply not truncating the Taylor series: 

<^^^nad2 = SP2 - clSp2 - -TT^Spi 



op 

r.2 



<5Pnad3 = SP3 - clSpz - i^Sp2Spi ~ ^Spi^ 



dp ^^"^^ dp^ 



(5.10) 

(5.11) 
(5.12) 



Thus, in the uniform density gauge, the pressure perturbation is equal to the non-adiabatic pressure perturbation at 
all orders. Then, Eq. (|5.9p becomes 



^ <5P„ad3 = 6(C2C( + C1C2) + 24C?C1 ^(C2'5/'„adl + C[SPn^2) ■ 



Po + Pq 



Po + Po 



(5.13) 



In the case of a vanishing non-adiabatic pressure perturbation, ^[ and C2 ^^re zero and hence we see that Cs is also 
conserved, on large scales. This was also found in Ref. [31, and previously in Ref. (23 |. 
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B. Klein Gordon equation 



The energy momentum tensor for a canonical scalar field minimally coupled to gravity is easily obtained by treating 
the scalar field as a perfect fluid with energy-momentum tensor [25| 



1 



Q/3 



where the scalar field is split to third order as 



and the potential U similarly as 



where we define 

SUi ^ U,,^5ipi 



SU2 = U^^^Sifl + U,^5ip2 



SU3 = U^^^^Sifil + 2U,^^5ipi5ip2 + U,^5ipz 



and making use of the shorthand notation [/_^ = Then, Eq. ()5.ip gives the Klein Gordon equation 



^0 V 



^0 



+ 8(f,^'„{H - 2H(t> - 2(f>') + B\ip'„B[ + 2Vi,» + + 2if'^^B\, + 0' - C'\{\ - 2</)) 

- 2C'^{B,,, - q,) + B'C^,,, - 2C'\,B, - B'<i>,, - 2B 
+ dipi,,{B'' + 2HB' + 2C'\j 

n.. . 

^' • 2nB' -2BjC'^' 



B' 



0} - 2C'^Sipi^,, 

+ 2nB') 

B'C^'j + B^Bj.' - B^B\j 



- B'B\, - 2C'^B', - 2C'W^k,j + '^C'^C.jk!" 



2C'^,j - 2B'> - - 2(j);'(f) - AHBjC'^ 



1 



,3 — 

+ 6dip2,rjC'^ - Q5lpi,^j{AC^'^ Ck - B'B') - 2An^'a(l){l - 20 + 40^ - 0' + 30'. 



- 125' ( Vi,. + 2-5^2 J + 24Vi,,(i3' 



6<po (20(1 - 20 + 402) + - 4B'B,0 - 2B'BJCy) - 6&\jip'o{Bi(j) - 25;) 



+ 6C"^^[,(l-20 + . 



4Cy0 - 4CkjC\) - 12^'oC'' 2B,Qfe,^ + 2B^Cfe,j - S,- , - B.B'j + B,(f>,j - BkQ,, 



+ 2CikB 



2B, 



AHB,Bj 



6^0 



B' 



B' B', - B'B,B^ 



2nB'B, - 2B*B,0' - 8HB'B,4: + AB'cj)^, - 2B' 



AB' 



B' 



B'BjB\ 



AB'B' 



0. 



(5.14) 

(5.15) 

(5.16) 
(5.17) 



(5.18) 



One can again see the coupling between first order tensor and scalar perturbations. For example, the SipijC^^ (jfj 
contains a term that looks like i5</3i,i/ii*-'0i j , which occurs only at third order and beyond. 

Again, we refrain from splitting up the perturbations order by order for ease of presentation. Once split up, one 
can then replace the metric perturbations by using the appropriate order field equations. We present the Einstein 
tensor at third order in the next section. Note also that Eq. (|5.18p implicitly contains the Klein Gordon equations 
at first and second order. We refer the reader to, for example, Rcf. for a detailed exposition of the second order 
Klein Gordon equation. 
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C. Einstein tensor 



The Einstein tensor, which describes the geometry of the universe, is defined as 



(5.19) 



where R^^^ is the Ricci curvature tensor and R is the Ricci scalar. Here, we give the components of the Einstein tensor 
up to third order: 



,1 



,1 



+ 2n{C',^ - B,,,)] + C,kAC'\' - ^C^'^;) + C\{B,,i - , + ^m) + 2G*^'.,G,fe. 
+ ^Bj4B\i + B\') - m'^m'' - B,B') - ^B\,Bj,^ - mB\,(j) + G% , 



■^'^jk,i L' k,ij ^ V '^i- 
k 



(5.20) 



+ B^ {Ckj^i' — C^kAj + Cik^k'' — ^"^Cij — 2'HBj^i) 



• 2G'^^' 
1 



G 



ik,j ~ ^jkA + 2 i^k,!] - Bi,kj) 



Bj/ - V^B, + 2 (g;./' - G^',, ,)] - C'^'Ckj,, + , 



(5.21) 



a^ff, = G^% + 27^G^ - \{B^\, + %*') - C\/ + G^,/ - 0/, - V^G^, + Q^/' - H {B\, + S,/) 

+ S', I (^H' - (1 - 20) + 2H (B\k - C^'k + 0') + B'^'k - C^'m - C^'^ + (0 + G'oj 



C:(fc/ + fcj — 2G' j_fe + 2T-i{Cjk^ + G\.j — G*j_a:) + - (-Bj/fc + B\jk — 2Bk.^ j) 



+ (G*^; - - B^fc)(G';■ - i {b\, + b,:)) + g"^' {B,.k - 2c[^) + g;,.b',^- + 0,^,, 



+ (B'^' - 2G''^i + GV*'- + 0,'-)(Gjfe.'' + C\^j - G%- ,) + -B\Bk/ - V^B, + 4^0,, - 2G"'fe,, + 2G;^>-) 



+ 2G*'' 
1 



{b'^m + 4j) 0,,fc - Gfc,,/ - C,iJ + V^Gfc, + G^.,fc + n (B,.k + i?.,, - 2G;,) 



{Bk:B\, + B,, «B\fc + , 



'~ik I /^l k/ 



(B^; + B\j + 20^ + 2H{Bj,' + S\j-) - 2G'j - AHC ^ 
+ (5', I (n^ - (402 - BkB^") + 20[G^"fc - B'^'fc - V^^ + 27^(G^'fc - 20 - B^,-) 



B' 



2C' i,k - 2G;,/ + V^Bk - Bi,k' + 2n{B'k - 0^^ - 2C'kj + G^fe) 
Gj.; — V^Cki + 2HCi^i + 2Cir,i,k"^ — C"^m,ki — 2HBi,k — Bij. — 



2B''' {C\.k ~ Cki!) 



C'''k[B\i--^ I 



1 



k-B\k] --[2B\kBu 



; I + .k^lm. + LyZm.fc I . " -L^ . I " L' I m. ^ -L' 



- BukB"".' - 3B\kBi.'') + 0,fc {C'l'' - 2G'^^ - 0,'-^ + G'dj^ , 



(5.22) 
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where G^o,G^i,G'^j are the third order corrections (the latter spHt into a diagonal part G^j, and an off diagonal part 
which we give in the appendix as Eqs. (|Aip . (|A2p . (|A3p and (|A4p . respectively. Note that, in calculating the 
third order components given above, we have implicitly obtained the full second order Einstein tensor components 
for fully general perturbations (i.e. including all scalar, vector and tensor perturbations). We shall extend this second 
order analysis, by presenting all the geometric and matter tensors as well as conservation and constraint equations in 
full generality, in a future publication pvj . 

VI. DISCUSSION AND CONCLUSIONS 

In this paper we have developed the essential tools for cosmological perturbation theory at third order. Starting 
with the definition of the active gauge transformation we have extended the work of Ref. [ll| to third order, and 
derived gauge invariant variables, namely the curvature perturbation on uniform density hypersurfaces, (3, and the 
density perturbation on uniform curvature hypersurfaces. We also relate the curvature perturbation ^3, obtained 
using the spatial metric split of Ref. [l3| to that introduced by Salopek and Bond [22|| , which is also popular at higher 
order. 

We have then calculated the energy and momentum conservation equations for a general perfect fluid at third order, 
including all scalar, vector and tensor perturbations. The Klein Gordon equation for a canonical scalar field minimally 
coupled to gravity is also presented. We highlight the coupling in these conservation equations between scalar and 
tensor perturbations which only occurs at third order and above. Finally, we have presented the Einstein tensor 
components to third order. No large scale approximation is employed for the tensor components or the conservation 
equations. All equations are given without specifying a particular gauge, and can therefore immediately be rewritten 
in whatever choice of gauge is desired. However, as examples to illustrate possible gauge choices, we give the energy 
conservation equation on large scales (and only allowing for scalar perturbations) in the flat and the uniform density 
gauge. This gives an evolution equation for the curvature perturbation ^3, Eg. (j5.13p . As might be expected from 
fully non- linear calculations psj and second order perturbative calculations [29j . the curvature perturbation is also 
conserved at third order on large scales in the adiabatic case. It is worth noting that higher order perturbation theory, 
as discussed in this paper, has the advantage of being valid on all scales whereas fully non-linear methods, such as 
the 5N formalism, are only valid in the large scale limit. 

Another application of our third order variables and equations, in particular the Klein Gordon equation (|5.18p . is 
the calculation of the trispectrum by means of the field equations. Whereas calculations of the trispectrum so far 
derive the trispectrum from the fourth order action, it should also be possible to use the third order field equations 
instead. The equivalence of the two approaches for calculating the bispectrum, using the third order action or the 
second order field equations, has been shown in Ref. [s^. Having included tensor as well as scalar perturbations it 
will be in particular interesting to see and be an important consistency check for the theory whether we arrive at the 
same result as Ref. [3T| . 

A final advantage of extending perturbation theory to third order is that, in doing so, one obtains a deeper insight 
into the second order theory. Also second order perturbation theory, despite remaining challenging, becomes less 
daunting having explored some of the third order theory. 
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APPENDIX A: EINSTEIN TENSOR 

Here, we give the third order corrections to the components of the Einstein tensor. We do not split up perturbations 
order by order. 
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,1 



2B^{Cki,i ~ ^ji,k) 



+ B,,,){4>B^: + BV/) - 2'HE'B^B,^, + 



BX2afc,'= - C\-,, + 0,, - 2HB,) + B\,0 - 2G\ 
Bk,jB\&\, - 2C,^^') + 2HB'' + 20(2Cy/' - C-''^- , + 20,, - 37^5,)] + 

B*B^(2Qfc,> - C^,,, - V^Cy ) + - 20,^!^ ~ 0,,) + B'^{B\,, + 2C^/' - 2C^"; ^ - V^B,) , (Al) 



— Bij){C''i^k — Cki!") + (Bi^i — Bij){Cjk,^ — 2Ck ,j) + 2{Bjj — Bi^j){C\k^ — Ci\k) + (Bj^i + Bij)4>^k 



+ 2{C'^k - Bk.M. + 20(2C^,,, - 2C[k^^ + - i?^,,,) 



2B, 



C^-' {C^i^ik — Cki,i — Cii^k + V'^Cik) 



+ c^'iCki/ - c^l,^k - aLk' + c^m) + ic'i," - + c^k,' - c^,^k) + CuAc'!' - 

+ \CkuC'''^ + 2C^0(C\,J' - 2C'=-'',fc) + B^' [B,{C[k.'' ~ C\, + ^{Bk/ - V^B,)) + CikB,' 



1 



+ B,4B\k - C'k) + 2Sj,fc(2Cr - 5^,5,,") 



20 



(S,,, - B,,){C''\k - ^C^k^^) + 0,J'(B,,, + B,,,) 



+ C^,C'=-'',, + 2C-''', 



2CU^ -BkAC: 



J k 



k j^ 



2i?,,^ 



+ 2CfcjC^''(i?,,,'=-i?,,,'=) 



0,,(C7\,, - 2Qfc,'= + + 20(C^,,, - CJk^. ~ C,k/ + V^Cy + WB,,,;) + 2H(B'=Qfe,, + 2Cfc,.B'=,,) 



+ 0,^(a,,fc - C,k,3 + C^jfc.*) + 20'(2C^';. , - 2C^,^■ - B,,,^' + V^B, - 12H0,,) ^ 



(A2) 



rt/^l/f jzyll Tit ll r)y^ra I 



2C' 



2Cj„,'™ - 2V2c'j - 20,/ 



■^jm, ^ rn., j \ ^ ra.j ^ j., mj 



+ 4C* C ''^\Ckm,jl — Clrn,kj + Cjm.kl — Cjkdm) + 4C ^'^CkliCj. 

- 2C7;.,(C'; - B\i - 0') + 2C7;,(2C'^ - Bk!) + C'[{Bk,, + B,A - 2C^,B,,' - (0' + + 



+ Bi,jB\k + S,, ,Bfc,' - 20,fc0,, - 20(B[,^- + B'^^k + ^.jk + 4HCfcK2C'; - B',^ - Bj!) - AHBk<j>„ 



C 



^iL- m,; — ^L-im, jlL- j,fc — O k,j — ^jk, ) + Z(ZLy — L-kl, )['^jm, + L- m,j — L- jjn) — iOkm, ^ l,j 

+ HCji^m — C'im,;)(C,fc"' — C^kJ"') ~ 2C*jCfc; + iC^'kCji + + Bj,^){Cki — Bi_k) — 20,'(Cj/c.* + C'^.j — C'^j,k) 



ZL- jOi^k — ZL^ k^jd ^ '^'-'jk^ ,i + ^k, — ^^k\'-'jl, + ^ Lj — j.l 



- B^Bk 



B'B^' 



B^B^ 



J - n{B\, + B,: - 2a'; 

f^'kL j — ^jl, fe — O IJk + O j^kl 



C\jk — C'kji — C'j,fei + V'^Cjk + 2HBk,^ 



B' 



c^kic- 



■I k 



20" 



CkiC^^j + ^eiBk^.+B^.k) 
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1 

2 



+ C\{AC"^j - 2C% - C,//' + 2B/,/ - Bjf - + C*;(2Cm; -B'f,~ C\.k + 0,fc) + 2C*"(Cj7,fe - Cj^^i) 
+ 2Cj,(C^',fc - C\}) + + C',,fc - 2Cfc,/) (B,; + B\,) + (C',/ + - C\j}){Bk.i - Bim - 2C'^i - mCki) 



20 



+ C"fc(2C'=;- - Sj- '=) + 2C*'^'Cj'fe - Cj;^,^','^ + -{B\j + Bj:)[C^'f, - 20' - B^fe) + C"j{B\k - C^'^ + 20') 



(27^5* - 0/) - 4Ha''=(2C;^. + Bk,j + Bj,k) + ^ {Bk^B^, + B\kB,,'') 



(A3) 



B' 



+ c 



rfci 



40,„„(^0 ^kl — ^^l ,k ) + ^^km\^ ,1 + ~ oO i) + Z(Ly „_ — /lO ^n)(^'~^lm,k — ^kl,m) 

+ C'A;„,m(6C/",'" - 2C;",") + 4C;m,fc0,™ - 2C""m,fe0,i - 2Cfei,m0,™ + 20,fc0,i + 4C/cm,™(0,i - C;„,") 

+ 2(i?,.fc - C^,)(S",™ + 0' + 47^0) + 40(0,,, + B,',, - C'li) + i?,„,fe(i?, ™ - 
+ 2Bk{2Cir„,'' - 2C":„_, + B"\un - V'-Bi - 2'H{B[ + C"'m,i - 2Q™," - 0,,)) - ^Sfe,™Bz ™ 

2C"'(C,'„, - Bi^,n) + cj)C"i - 8H02 - 400' - 20S',i + B\b'i + C"\,n,i - 2C/™," - 0, + 2HB,) 



+ B'' 



Bk ( C';' - B"i - V^4> - 2W(20' + B\i)) + (2G,„,fc'" ~ C™™,,; - V^C^ - 2H(C^, + B/,,) 



2a" 

a 



I (2Cfc; + 40) + 2C''"( 2C[.„ ; — 2C;'„ J, + Brn,lk — -Bfc,mi — 21-L{Clm,k — 2Ckrn,l 



c 



•kii 



- B^'^fc - V20 + 2rL{n4> - 30' - B\k)\ + C^"0(2B,,,- - 3C^, j + -Bi,k4>{B\' - 3B',^) 



+ B^fc0(40' + B\i) + 20(Cm, V/ - CVfc0,' + 20,fc0,'=) + 2B'='(2C,,/0 - 0C',,fc) 



(-Sfc,i — Si,fc — 2C^;) + C'™,m(4C^; — 2Bkj) + B\i{2Ckm™ - C™'m,k - B'k + 0,fe) 



+ B™,z(C'fc;" - 3C"fe.' + C"",fc) + C[„{Ck"\' - C"",fe) - 0'(C'z,fc - 2Cm/) + i3',fcani, 
+ 20(2Cfci/' - 2C% + BiJ - V^Bk - 4H(2S^ + C^fc ~ 2^,;/ - 20,,) 



(A4) 
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